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Abstract 

As a continuation to our previous work [9j [10], we consider the 
domino tiling problem with impurities. (1) if we have more than two 
impurities on the boundary, we can compute the number of corre- 
sponding perfect matchings by using the hitting matrix method [4]. 
(2) we have an alternative proof of the main result in [9] and result 
in (1) above using the formula by Kenyon- Wilson (6j [7j of counting 
the number of groves on the circular planar graph. (3) we study the 
behavior of the probability of finding the impurity at a given site when 
the size of the graph tends to infinity, as well as the scaling limit of 
those. 

Mathematics Subject Classification (2000): 82B20, 05C70 

1 Introduction 

1.1 Background 

Let G = (V(G), E(G)) be a graph. A subset M of E(G) is called a perfect 
matching(or a dimer covering) on G if and only if for any x G V(G) 
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Figure 1.2: An example of Gi,t for k = 2. 

there exists e G M uniquely with x G e. Let .M(G') be the set of perfect 
matchings on G. The perfect matching problem was first studied by Kaste- 
leyn, Temperley-Fisher, in the context of the statistical mechanics, and many 
papers appeared since then. In particular, for the perfect matching on the 
bipartite graphs, such as domino or lozenge tilings, many results have been 
known (e.g., [5j and references therein). 

In this paper we consider the perfect matching on a non-biparite graph 
defined below. We use the same notation in |TU]. Let G 2 be a subgraph 
of the square lattice, and let G\ be its dual graph(Figure II. ip . 

Let G\ t T be the graph made by adding (2k — 1) vertices T 1; T 2 , • • • , T 2 k-i, 
which we call the terminal, to the boundary dG\ of G\ (Figure fL2l) . 

We next superimpose G 2 and Gi,r, and we make a vertex wherever two 
edges cross, and call it a middle vertex. For the vertices in the boundary 
of G\, we add extra edges toward the outer face of G 2 , and make vertices, 
which we call the boundary vertex, wherever two edges meet (Figure II. 3p . 

By adding furthermore diagonal edges alternately, we have the graph 
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Figure 1.3: We superimpose G2 and G\ : t and then add middle and boundary 
vertices 
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Figure 1.4: Graph G (2) 



G^) (Figure OD- 

Figure [L5l shows an example of a perfect matching on G^ 2 **, where we have 
two kinds of edges. For M G Ai(G), we say an edge e 6 M is an impurity if 
it is an edge between the vertices of G^t an d Gi- The number of impurities 
is constant on M.(G^) and is equal to k. 

We shall review the known results on the perfect matching problem on 
G^ k \ In [8], they consider two types of elementary moves which transforms 
a perfect matching on G^ k ' to a different one, and showed that any two 
perfect matchings are connected via a sequence of elementary moves. Then 
we can construct a Markov chain on Ai(G^) with the uniform stationary 
distribution. By MCMC simulation, we conjecture that all impurities tend to 
distribute near the terminals, more precisely, the number of perfect matchings 
is maximized when all the impurities are on the terminals. Ciucu [TJ |2] 
studied the dimer-monomer problem mainly on the hexagonal lattice and 
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Figure 1.5: A Dimer Covering on G^ 

found that monomers interact as if they were the charged particles in 2- 
dimensional electrostatics. Though his setting of the problem is different 
from ours, his works should have something to do with our conjecture. In 
[H], they studied the one impurity case, and derived a formula for computing 
the number of perfect matchings if the impurity is arbitrary fixed. From 
that formula, the conjecture above is generically correct for k = 1, in the 
sense that the probability of finding the impurity on the fixed site (under 
the uniform distribution on Ai(G^)) decays exponentially away from the 
terminal. 

In this paper, we study the k > 2 case and derive a formula to compute 
the number of perfect matchings if the impurities are fixed on the bound- 
ary (Theorem ll.2p . Moreover, we study the behavior of the probability of 
finding the impurity when the size of the graph tends to infinity. 

Notation 

We collect the notations frequently used in this paper. 



(1) A G is the Laplacian on a weighted graph G : 

(A G f)(v)= £ c vw (f(v)-f(w)) 



where c vw is the weight on the edge (v, w) G E(G) which we take c vw = 1 in 
this paper. 

K is the Laplacian on Gi t restricted on G\ : 



K :— A GlT \ Gl — l Gl A GlT l Gl . 



;i-2) 
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Figure 1.6: G^t for one impurity case 



Is is the characteristic function on a set S. 

A(x,y) or A x y is the matrix element of a matrix A. 

(2) jW, l( 2 \ ■ ■ ■ , G E(G^) are the location of impurities, where I® = 
4 j) G V(Gi), j = 1,2,- --,k, i = 1,2. For the one impurity 
caseO = 1), we simply write I = (I U I 2 ) e h G i = 1,2. 

ti,t 2 , ■ ■ ■ ,t 2 k-i G ^(Gi) are vertices of G\ connected to the terminals 
Ti, T 2 , ■ • • , T 2 fc_i respectively(Figure II. 2p . For one impurity case, we simply 
write t (Figure II .6p . 

In the following subsections, we summarize the results obtained in this 
paper. 

1.2 One impurity case 

We first recall the result in [H] where we set k — 1 and counted the number 
of matchings when the location of the impurity is given. 

Theorem 1.1 Let x G G\. Then the number M(x) of perfect matchings on 
G^ whose impurity satisfies I\ = x is equal to 

M{x) = {(K-^fat) detK\. 

If we specify I± = x G G±, we have four possibilities of putting I 2 G G 2 , 
but the number of perfect matching is independent of the choice of that. In 
Section 3.1, we give an alternative proof of Theorem II .11 by using the theory 
developped by Kenyon- Wilson [61 [7] . 
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Figure 1.7: with fixed impurities 

1.3 Two impurities case 

Set k = 2 and let a, b G Gi be vertices on the boundary dG\ of Gi such 
that there are no terminals in between (Figure [L~Tj) . We put the impurities 
such that their location jW,/^ satisfies /j = a, ij = 6 and the other 
ends I2 j ^2 ^ e 011 the boundary <9G2 of G2. We remark that, unlike the 
one impurity case, the number of perfect matchings becomes different when 
we put /W,/( 2 ) ^ <9G 2 (Figure [L8|) . Let dC be the set of vertices on the 
boundary dGi of G\ which lie between a and b (Figure [L9|) . 

Theorem 1.2 Suppose the two impurities satisfy l[ = a, 1^ = b and 
Jjj G j = 1,2. Then the number M(a,b) of the corresponding perfect 
matchings is given by 







f 


L a ,ti 


La,t2 


La,t 3 \ 


M(a,b) = A(a,b) : = 


det 






Ldc,t 2 


L 9C , U \det{K) 






{ 


Lb,tx 


Lb,t 2 


L b,t 3 ) 



where 

L xM = (K~ l )(x,ti), x = a,b, i = 1,2,3 

L dCM = ^(K-^U), i = 1,2,3. 

yedc 

This determinantal expression implies that impurities are repulsive each 
other. We give two proofs of Theorem 11.21 : one is due to the hitting matrix 



6 





-o- 


--Q- 


--Q- 

i 




-Q- 




-<> 


6- 


-6- 


-6- 


--6- 

i 
i 


M 


-6- 


-6-- 


i 

-o 


6- 


-6- 


-6- 


--6- 



Figure 1.8: This impurity configuration is not allowed for Theorem 1.2. 




Figure 1.9: G\ t for two impurities 
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Figure 1.10: The case where Gi is an one-dimensional chain and the corre- 
sponding G^ . 

method by FominJI] in Section 2, and the other one by Kenyon- Wilson [61 [7] 
in Section 3.2. If 1% ^ (as in Figure [L8]) . then we have another for- 
mula given in Section 3.4. Essentially the same formula also holds for k > 3 
case to be shown in Section 3.5. However, if some impurities are not on the 
boundary (i.e., 7^ ^ dG\ and 7i; ^ dG2), then our methods of proof do not 
apply and it would be difficult to compute the number of perfect matchings. 
On the other hand, if G\ is the subgraph of the one- dimensional chain(Figure 
ll.lOp . we can compute the number of perfect matchings for arbitrary position 
of impurities(Section 3.6). 

1.4 Large size limit 

In the one impurity case, we consider the limit of M(x) as the size of G\ 
tends to infinity, for the following two graphs. 

f Gf{n) := {(x,y) \ x = 1,2, •• -,n, y = 1,2, ■ ■ -,n} (2-dim) 
1 \G?(n) := {1,2, ■■•,n} (1-dim) 

We can also study the case where G\ is the n x m-rectangle, provided 
n, m grow proportionally. We set the problem for the 2-dim case below, 
but the 1-dim case is formulated similarly. Connect the terminal T to 
r = (1, 1) e Gf (n) (for 1-dim, r = 1 G Gl h (n)). We consider the following 
two problems. 

(1) Compute lim^oo P (7i = (x,y)) for fixed (x, y) G G{ q (n). 

(2) We scale G S i(n) by - so that it is contained by the unit square, and 
consider the probability of finding the impurity on a small region in it. In 
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other words, for any fixed < c\ < c 2 < 1, < di < d 2 < 1 we would like to 
compute 

/i([ci,c 2 ] x [di,d 2 ]) := lim P(/i G [cxn,c 2 n\ x [c^n, d 2 n]) 
of a measure /i on [0, l] 2 . 

For the first problem, 

Theorem 1.3 (1) (Example 3.6 in [9]) For the one- dimensional chain 
Gl h (n), for fixed j G Gl h (n), we have 

P(h=j) = \\+ j (l + o(l)), n-^oo 

where X + := 2 + a/3. 

(2) For the two-dimensional grid Gl q (n), for any fixed (x,y) G G'l 9 (n), we 
have 

limP(/ 1 = (x, 2 /))=0. 

Theorem [L3]has the following implications, (i) for the one-dimensional chain, 
the probability of finding the impurity decays exponentially and it is localized 
near the terminal, (ii) for the two-dimensional grid, the probability spread if 
the size of the grid is large. 

The perfect matching on G^ is determined by a spanning tree T on a graph 
Gi } r defined in Section 2.1 (Theorem 12. lj ). And the difference between the 
one- and two-dimensional cases comes from that of the expectation value of 
the length of T (Proposition 14. ip . In one-dimensional case, it is bounded 
with respect to n, while it diverges in the logarithmic order in the two- 
dimensional case (Proposition 14.1( 2)). This observation together with the 
Chebyshev's inequality also solves the second problem : 

Theorem 1.4 In both cases, fi is equal to the delta measure on the origin. 

In [3] , they studied the lozenge tiling with a gap and showed that the corre- 
lation function behaves like =, where r is the distance between the gap and 
the boundary. In our case, the impurity tends to be attracted to the terminal 
T and not to the whole boundary, so that the situation is different. In the 
following sections, we prove those theorems mentioned above. 
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Figure 2.1: G^t,r '■ all vertices on dG^T are connected to R. 

2 Hitting Matrix 

2.1 An extension of Temperley Bijection 

We first recall the results in \X0\ . We consider an imaginary vertex R (called 
the root) in the outer face of G^t and let G\,t,r be the graph obtained by 
connecting all vertices in dG\pp to the root -RfFigure 12. ip . 

Tl-tree is a tree on G^r starting from the root R, which directly con- 
nects it to a terminal, and ends at a vertex of G\ t T- TO-tree is a tree on 
Gi,t,r starting at a terminal and ends at the root R through a boundary 
vertex. IO-tree is a tree on G^t,r starting at a vertex of G\ and ends at 
the root R through a boundary vertex(Figure I2T21) . Theorem 2.8 in [10], in 
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TO-tree 




Figure 2.2: Example of Tl-tree, TO-tree, and IO-tree. R is omitted, 
a slightly different form, is : 

Theorem 2.1 We have a bijection between the following two sets. 

M(G^) := { perfect matchings on G^} 

Q) := {(T, 5*, {ej}j =1 ) | T : spanning tree on G^t^r, 
S : spanning forest on G2, 

{ej}j =1 : configuration of impurites, with condition (Q) } 

Q 

(1) T is composed of k Tl-trees, (k — 1) TO-trees, and the other ones are 
10 -trees, 

(2) S is composed of k trees, 

(3) T , S are disjoint of each other, and the k Tl-trees of T and the k trees 
of S are paired by impurities. 

Figure 12.31 is the spanning forest S on G2 corresponding to the spanning 
tree of G\^t,r shown in Figure 12. 2[ and Figure 12.41 shows the corresponding 
perfect matching of G^ . Sometimes Tl-tree and TO-tree sticks together to 
form another tree(Figure 12. 5j) . For one impurity case(/c = 1), T is composed 
of a Tl-tree and some IO-trees and the bijection in Theorem 12. II is the same 
as Temperley's bijection(Figure IT76]) . 
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Figure 2.3: The spanning forest on G2 corresponding to the spanning tree 
in Figure 2.2. 
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Figure 2.4: The corresponding perfect matching on 

TI + TO-tree 




Figure 2.5: TI + TO-tree 
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Figure 2.6: The case of one impurity 

2.2 Proof of Theorem 11.21 using the hitting matrix 
method 

By Theorem 12.11 M(a,b) is equal to the number of elements in FiG^^Q) 
such that the corresponding spanning tree of G\,t,r is composed of (i) a 
Tl-tree connecting a and Ti, (ii) a Tl-tree connecting b and T3, and (hi) a 
TO-tree connecting an element of dC and T 2 : 

M(a, b) = I F is a spanning tree on G\ } t,r composed of 

(i) Tl-tree connecting a and Ti, (ii) Tl-tree connecting b and T 3 
(hi) TO-tree connecting dC and T 2 , and 

(iv) other IO-trees } (2.1) 

Let Gi t R be the dual graph of G 2 in the sense that we regard the outer face 
as a face of G 2 and put a dual vertex R there, which is called the root (Figure 
12.71) . Gi.r is also obtained by identifying all terminals in G\ y T,R with R. 

Then the RHS of ( 12. ip is equal to the number of spanning trees T of G\^r 
with the following property : cutting all the edges connecting to R, yields a 
spanning forest Ti of G\ which satisfies (i) a and t\ are connected, (ii) b and 
£3 are connected, (hi) dC and ti are connected : 

M(a, b) = {T I Ti is a spanning forest on G\ such that 

(i) a is connected to t\, (ii) b is connected to £3, 
(hi) dC is connected to ti, } 

By taking acount of Wilson's algorithm [TT] generating the spanning trees 
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I [ 

Figure 2.7: In G\ t R, all vertices on the boundary of G\ are connected to R 
such that their degree = 4. 



on G\ % r uniformly at random, we have 

M(a,b) 



(j{ spanning trees on G^r} 

= J2( s § n °) w{-K{a,t ai ))w(Ti(b,t (J2 ))w('K(dC,t a , A ))l A{a) . 
(j 

In the notation above, tt(x, y) is the walk from x to y, |tt| is the length of the 
walk 77, w(ir) := 4~' 7r ' is the weight of the walk it, and 1a(<j) is the indicator 
function of the conditition A(a) : " Tc(b,t a2 ) is disjoint from the loop erased 
part of 7r(a, £ CT1 ) and 7t(c, t^) is disjoint from the loop erased part of tt(6, £<t 2 )" 



A((t) := {<ir{b,t a2 )nLE(ir(a,t ai )) 



7T 



( C ,t CT3 )nL£(7r(Ma 2 )) =0} 



Therefore by Theorem 7.2 in |4| and the matrix tree theorem : 
jj{ spanning trees on G^r} = det(K) we have 



M(a,b) 



( H. 



det 



a,ti 



a,t2 



b,t 2 



H, 



a,tz 



det(K) 



\ H 9CM H d c t t 2 Hq C . 



t 3 



(2.2) 



where H x ti is the hitting probability from x to t\ : probability of a simple 
random walk on starting at x hitting the root R through t iy and 



H. 



oc,u 



zedc 
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To compute H x y , let Q be a matrix on G\ given by 



:= [\ ( x ~ v) 

\ (otherwise) 



Q is the transition probability matrix of the SRW on G^r restricted to G\. 
Since 4(1 — Q) = K, we have 



for x G Gi, i = 1, 2, 3. Substituting it into (Q yields Theorem O 

3 Alternative proofs by a Kenyon- Wilson's 
result 

3.1 Kenyon- Wilson's results 

We first recall the results by Kenyo n- Wilson [6j E]. Let Gc be a planar, 
finite, and weighted graph, and for e = {v,w} G E(Gc), let w(e) = c vw be 
the weight on that. We choose n vertices on the boundary of Gc, number 
those points anticlockwise, and let iV := {1,2, ■ ■ ■ ,n} which is called the 
node. A graph with such a property is called a circular planar graph. 
Let / := V{Gc) \ N. We say a spanning forest g on Gc is a grove if each 
tree composing g contains at least an element of N. A grove g determines a 
planar partition (non-crossing partition) a of N. We set 



For a grove g, let w(g) := Yie^g w(e) be its weight. For given planar partition 
a G V p i of N, we consider all groves which realize the partition o and let 
Zc c (cr) '■= Sgeo-MsO be their weighted sum. 

According to the decomposition V(Gc) = N U I, the Laplacian Ac c on 
Gc can be written as 



H x ,u = ((I - Q) ^)x,u ■ j 
= K Xjt . = L x>ti 



V := { partitions of {1, 2, • • • , n}} 
V p i := { planar-partitions of {1, 2, • • ■ 



n}} 




(3.1) 
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We define the response matrix L by 



L = GK H - F. 



L has the following physical meaning. We regard Gq as an electrical circuit 
with c vw being the conductance on the edge {v,w}. For i,j G N, set the 
voltage / on iV by 



Then L it j is equal to the current flowing out of j. 

a G V p i is said to be bipartite if it satisfies the following condition : 
(i) it is possible to color the vertices in N contiguously into red and blue 
such that only the nodes of different colors are connected, (ii) if a node lies 
between red region and blue region, it can be split into two colors, and (hi) 
the number of each part in a is at most two. 

For bipartite a, let Lx,y be the submatrix of L such that the color of the 
row(resp. column) is X (resp. Y). We then have the following theorem. 

Theorem 3.1 (Theorem 3.1 in [7]) 

For bipartite a G V p i, we have 



Remark 3.2 Theorem 3.1 in /?]/ treats more general case called the tripartite 
partition and write Z(a) in terms of Pfaffian. Theorem \3.1\ also holds for 
non-planar graph G provided G satisfies the following condition : 

(A) For any non-planar partition r G V\V v i, we always have Zg{t) = 0. 

3.2 One impurity case alternative proof of Theorem 




(otherwise) ' 



x = i 



Z{a) = IdetL^B • det(i^)|. 



For x G Gi, we define a matrix K x by 




K(x,x)+1 ((y,z) = (x,x)) 
K(y, z) (otherwise) 
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Figure 3.1: Gq in one impurity case : add X to x, identify all boundary 
vertices with Y 



Lemma 3.3 



M(x) = \(K- 1 )(x,t) det(K x )\. 



Proof. We construct a circular planar graph Gq by the following proce- 
dure(Figure !3.ip : (i) add a new vertex X to x, (ii) add all boundary vertices 
to Gi, and identify them. Let Y be the corresponding vertex and we set 



Although Gc is not planar, the condition (A) is satisfied, for iV contains 
only three elements so that we do not have non-planar partitions. Let a = 
XT\Y be a partition such that X and Y are connected and Y is isolated. 
By Theorem 12. 1[ we have M(x) = Zq c {<j). Writing Ag c as in ( 13. ID . the 
response matrix L satisfies L = GK~ l H — F. Note that we added an extra 
vertex X to x, so that we have to use K x in (13. ip instead of K. Since F = 0, 



V{G C ):=IUN, I 



V(Gt), N:={T,X,Y}. 



L X , T = (GK~ 1 H)(X,T). 



Using 



(H) 



(G) 




-1 (i e N,j e I,i ~ j) 

(otherwise) 

-1 {i E I,j E N,i ~ j) 

(otherwise) 



17 



and the fact t is the only vertex connected to the terminal T, we have 

L x , T = (K x )- 1 (x,t). 
Substituting it to Theorem 13.11 yields the conclusion. [] 

Proof of Theorem \l.l\ We show that the dependence of Lx,t = {K^^^x.t) 
on x cancels with that of detK x . In fact, by taking the (x, x)-component of 
the resolvent equation 

K- 1 - K- 1 = K~\K - K x )K~ l (3.2) 

we have 

{K^){x,x) - (K~ 1 )(x 1 x) = -(K~ 1 )(x,x)(K~ 1 )(x,x) 

thus 

= IW& (3 ' 3) 

Moreover taking the (x, t)-component of (13. 2p and using (13. 3|) . 

Let k x be the submatrix of K x by eliminating the x-th row and x-th column. 
By Cramer's formula, 

( Tr\-H \ det(^) 
{K) (X ' X)= deIM' 
On the other hand, expanding the determinant yields 

det(fQ = det(K) + det(k x ) 

leading to 

It then suffices to substitute (j3,4fl . (13. 5p to Lemma [3~31 to finish the proof of 
Theorem 11.11 [] 
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Figure 3.2: Gc in two impurity case 



3.3 Two impurity case I alternative proof of Theorem 

As we did in Section 3.3, we construct a circular planar graph Gc as fol- 
lows(Figure 13. 2j) : (i) add all the boundary vertices to G±, (ii) Let A, B 
be the boundary vertices attached to a and b respectively, and identify all 
boundary vertices lying between A and B, which we call C, and (iii) iden- 
tify the other boundary vertices in their respective region. Let D, E, E, G be 
these vertices and set 



V{G C ) ■= I U N, I := V{G X ), N := {T h T 2 , T 3 , A, B, C, D, E, E, G}. 

Let ti, t 2 , t 3 , a, b e I be the vertices connected to the nodes Ti, T 2 , T 3 , A, B 
respectively. By (12. II) . M(a,b) is equal to the number of the groves on Gc 
such that A <B- Ti,<9C o T 2 , B o T 3 and D,E,F,G are isolated. The 
corresponding partition a 6 P p / is bipartite (Figure 13. 3j) so that by Theorem 
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Figure 3.3: The corresponding planar partition o. P, Q, R are red, A, B, C 
are blue, and D, E are split into red and blue. 



we have 



» 



det 



L 



A,T 2 



L A 



■ T, 



( I- A. 

Lc,Ti l c,t 2 L c ,t 3 ] • det (if) 
V Lb,Ti Lb,t 2 Lb,t 3 



(3.6) 



It suffices to compute these matrix elements of L. Under the same nota- 
tion in flUJ), we have L = -F + GK~ X E = GK~ l H. Hence 



L A , Tj = {GK- l H) A , T] = (K-\, tj , j = 1, 2, 3 
likewise for Lb,t-- Similarly 



L c , Ti = (GK-^on = E (K- 1 : 

zedc 



z.t. 



(3.7) 



(3.8) 



By substituting (13. 7p . (13. 8p into (13. 6p . we complete the proof of Theorem 11.21 
□ 

Remark 3.4 If an impurity lies inside G\, we generally have a non-planar 
partition r with Z{r) ^ so that the above argument is no longer valid. 



3.4 Two impurity case II : impurities lying "near" the 
boundary 

In G^ 2 \ we consider the case in which both impurities satisfy l[ G dGi, 
j = 1, 2, but 1% ^ 9G2 (Figure l3T4"l) . By Theorem 12 .1[ there are two types of 
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Figure 3.4: If I 2 t dG 2 , we have extra spanning trees. 



corresponding spanning trees : one is given in Figure [3T5| whose contribution 
has been computed in Theorem 11.21 and the other one in Figure 13 .61 As is 
explained in Figure ll.6[ the tree connecting 7\ and T 2 in Figure 13.61 can be 
regarded as a composition of a Tl-tree and a TO-tree. In this subsection we 
compute the contribution from the spanning trees in Figure 13.61 Let c G Gi 
be the vertex connected to I 2 through a diagonal edge(Figure 13. 6j) . 

Theorem 3.5 Let a,b G dG\. The number M(a,b) of perfect matchings 



satisfying I[ 



(i) 



= a, l[ = b and I 2 4- 9G 2 , I 2 G dG 2 is equal to 
M(a,b) = A(a,b) + B(a,b) 

B(a,b) :-- 



det 



a,t 2 


La,t 3 ) 




b,t 2 


Lb,t 3 


■ det(K) 


c,t 2 


Lc,t 3 J 





where A(a,b) is the one given in Theorem \1.2\ and 



L^-^iK-^y). 



It is also possible to deal with the case where I 2 ^ dG 2 j = 1,2, with some 
extra terms computed similarly. 
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Figure 3.5: A spanning trees contributing to A(a,b). IO-trees are omitted. 
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Figure 3.6: A spanning trees contributing to B(a,b). IO-trees are omitted. 
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Figure 3.7: Gc to compute B(a,b). 



Proof. We compute B(a, b) which is the number of the spanning trees on 
G\ t T,R corresponding to the ones in Figure I3TB1 As is done in subsection 3.3, 
we add all the boundary vertices to G\, and let A, B, C be the boundary 
vertices connected to a, b, c respectively. We identify all boundary vertices 
between C and B, which we call D (Figure 13. 7J . Similarly, we define vertices 
E,F,G,H. Let 

V{G ) :=IUN, 1 = V(G X ), N := {T x , T 2 , T 3 , A, B, C, D, E, F, G, H}. 

B(a, b) is equal to the number of groves such that T\ 4-> A, T 2 4-> C, T 3 o B, 
and the other nodes are isolated. Since the corresponding a G V p i is bipartite, 
we can apply Theorem 13.11 [] 

3.5 More impurities case 

It is possible to extend our argument in Section 3.3 to the case where the 
number of impurities satisfy k > 3, provided all impurities lie on the bound- 
ary. Let If\ • ■ ■ , if ') = (di, a 2 , • • • , a k ), ctj G dG u j = 1, 2, • • • , k be 
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Figure 3.8: An example impurity configuration in one-dimensional chain 



the position of impurities. Let dCj, j = 1, 2, ■ • • , k — 1 be the set of vertices 
in dG\ lying between a,j and a J+ i. 



Theorem 3.6 Let k > 2. //a// impurities lie on the boundary and sat- 
isfy I? = a j} a, e dG 1 , I { 2 j) G dG 2 , j = 1,2,- 
corresponding perfect matchings is equal to 



M(a 1 ,a 2} ■■■ ,a k ) 



( L 



det 



LdCi,ti 



J ai,t2 



Ldc±,t 2 



We also have the analogue of Theorem 13.51 



k, then the number of 
\ 



J a\.t 2 k-\ 



L dd,t 2k 



det (AT) 



3.6 One-dimension case 

When G\ is the one- dimensional chain, the impurities always lie on the 
boundary. Hence we can compute the number of perfect matchings for all 
configuration of impurities, by using the methods discussed in Sections 3.3, 
3.4. For instance, we consider the case described in Figure [3781 There are the 
two types spanning trees on G^t,r given in Figures 13.91 13.101 For each case 
we construct the circular planar graph Gc as we did in Sections 3.3, 3.4. The 
number of spanning trees in Figure 13.91 is equal to A(a, b) given in Section 
3.3. The number of spanning trees in Figure [3.101 is equal to C(a,b), which 
is the number of groves satisfying T\ -H- T 2 , T 3 <H- b, and the other nodes are 
isolated. 
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Figure 3.9: A spanning tree contributing to A(a, b) 




Figure 3.10: A spanning tree contributing to C(a, b) 

4 Large scale limit 



In this section we go back to the one impurity case and prove Theorems II. 3| 
11.41 We take G\ to be the n x m rectangle : 

G x = Gi n,m) := 1 1 < % < n, 1 < j < m}. 

The terminal T is attached to the vertex (1, 1). Let It be the length of the 
Tl-tree of the spanning tree T given in Theorem 12.11 and let 

A (n,m) := ( Agit | Gi )-1. 



By Theorem II. 1| 
(1) P(Ii 



A( n < m X(x,y), (1, 1)) 



4EJ = iEli4, m ((x',2/'),(M)) + 2 



(2)E[Z T ]= E E^ (n ' m) ((^Z/'),(l,l))- 

x' = l y' = l 



(4.1) 



(4.2) 



We would like to study the behavior of P(/i = (x,y)) and E[l T ] as n = m 
tends to infinity. 
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Proposition 4.1 Let n = m. 

(1) For fixed (x,y) G Gf' m) , 



\im A^ n \(x,y), (1,1)) = A(x,y) 

, . . 2 f f sin(#:r) sm(6y) sin#sin</> , , 

a ■=— — -^r — i—^dew. 

it 2 J j(o,tt) 2 2 — cos — cos 
(%) Under the uniform distribution on the spanning trees, 

E[h]= E A^ n \(x,y),(l,l)) = -\ogn(l + o(l)). 

(x,y)eGi 7r 

as n — )> oo. 

Proof. (1) The eigenvalues and the corresponding normalized eigenvectors of 

j{(n,m) are gj ven by 

a ( kit \ ( lit \ 

e\.\ := 4 — 2 cos — 2 cos 



+ 1/ \m + 1/ 

2 . / kit \ . / Ztt 
w kl (x,y) := sin -x sin -y 

^/(n+l)(m + l) U+l ' 1 ' 
(A; = 1, 2, • • • , n, I — 1, 2, • • • , to). Hence 

/) rami 

= W — 

(n + l)(m+l)i^^e fcI 

/ /C7T \ . / lit \ . ( kit \ . ( lit 

x sin x sin y sin sin 



x 



n + 1 J \m + 1 J \n + 1 J \m + 1 

^°°A(x,y). 



n=m^oc 



Here we note that the behavior of the integrand near the singularity (9, <fi) 
(0,0) is 

sin(#x) sm(6y) sin^sinr/) 9 ■ 6 . .„ . . ., . , ,, N 

4- 2c o r/-2co S / = Sm(to) ' Sm( *» + 0(191 + 
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so that it is removable. 

(2) The sum of A^ m \(x, y), (1, 1)) is equal to 

E^>((s,y),(l,l)) 

x,y 

4 j_ sm 2 (^) sm 2 

In what follows, we let m = n. Due to the monotonicity of the integrand 
near the origin, we have, for a positive constant C, 

-£A n , n ((x,y),(l,l)) 

x,y 

> 1 /"tt rir sin 2 # sin 2 d6d<p 

hh (1 -cos0)(l -cos0) ' 4-2cos#-2cos0 _ 

Since the integrand behaves 

fl 2 2 1 4 

| . £ ' ^ 2 + 2 ~ ^ + 02 

near the origin, we have as n — > oo 

1 r 71 " n sin 2 # sin 2 dOdcf) 
7T 2 " 7i Ji (1 - cos0)(l - cos0) ' 4-2cos0-2cos0 

l /-t 4 rf 

7r Ji r Jo 
2 

= -logn(l + o(l)). 

7T 

On the other hand, 

^ 1 /■*" sin 2 sin 2 d6W0 

2-, y>-^]i. (1 - C os0)(l -cos0) ' 4-2cos0-2cos0 



+ 7- 

(n 



[21+11 

1 JL J sm —r am n - , 

n+l rn+1 

+ 1) 2 U (4-2cos^-2cos^)(l-cos^)(l-cos^) 
+(n «-» m, k I) 
I + 11 + III. 



, i — I • 2 it • 2 (2Z-1W 

4 1 JL J sm — -f-r sm ^ — 
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The second term is estimated as 

4 ^ (l + cos^Xl + cos^-^ 

j j J> n+l > v m+1 - 

l) 2 V 4- 

n+l m+1 

4 „ (n + l) 2 



+ 4 _2cos^ T -2cos^-^ 

' n+l m+1 

and similarly for Therefore the upper and lower bounds have the same 
leading order yielding 

E[l T ] = ]T A n>n {(x, y), (1, 1)) = - logn(l + o(l)). 

77" 

2,2/ " 

□ 

Proof of Theorem \1.3\ 

(1) follows from [9], Example 3.6, and (2) follows from eq. fl4.ip and Proposi- 
tionSUl), (2). □ 



Proof of Theorem \1.4 



Let I\ = r G G\ be the position of the impurity and take any c > 0. By 
Theorem 12. II if |r| > cL then 1? > cL so that by Chebychev's inequality 

P (Irl > cn) < P (Z T > cn) < —E[l T ]. 

cn 

For one-dimension, E[/t] < oo by Theorem 11.3( 1) and for two-dimension 
E[/ r ] = (const.) log n(l + o(l)), n oo by Proposition 14.1( 2). Thus 
lim^ooPQrl > cn) = 0. [J 
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